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Abstract 

The natural parameterization of vacuum oscillations in three neutrino flavors is 
studied. Compact and exact relations of its three parameters with the ordinary three 
mixing angles and CP violating phase are obtained. Its usefulness is illustrated by 
considering various applications: the study of the flavor ratio and of its uncertainties, 
the comparison of expectations and observations in the flavor triangle, the intensity of 
the signal due to Glashow resonance. The results in the literature are easily reproduced 
and in particular the recently obtained agreement of the observations of IceCube with 
the hypothesis of cosmic neutrino oscillations is confirmed. It is argued that a Gaus¬ 
sian treatment of the errors appropriately describes the effects of the uncertainties on 
neutrino oscillation parameters. 


1 Introduction 

After IceCube results, see e.g. [DEI El, the importance of a precise description of os¬ 
cillations has increased greatly. In the present paper we discuss a natural, easy-to-use 
and completely general choice of the relevant parameters. We illustrate its usefulness by 
quantifying the impact of the uncertainties on various physical quantities, implied by the 
imprecise knowledge on oscillations. 

We begin by recalling the main achievements in the discussion of cosmic neutrino oscilla¬ 
tions. The general formula for the vacuum averaged oscillations was given in JTj. Ref. [5] 
studied for the first time the implications of the observed oscillation phenomena on cosmic 
neutrinos. Various authors remarked the possibility to measure flavor ratios, possibly aim¬ 
ing to constrain the parameters of oscillations, e.g. mmmm- The relevance of oscillations 
for the interpretation of Glashow resonance was noted in |1U| . In [11] the single parameter 
that rules cosmic neutrino oscillations and depends linearly upon unknown quantities was 
identified; then, it was remarked |12j that this leads to a strong correlation between the 
effect of the oscillation parameters on the probabilities of oscillation, lessening chances of 
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measuring the oscillation parameters. Non-linear effects were studied in various subsequent 
papers including mmm- An interesting expression for all relevant parameters of vac¬ 
uum oscillations was proposed in m within a scheme of approximation aimed to improve 
the tribimaximal mixing matrix. In m m is] an expansion in second order of small 
parameters of this parameterization was obtained and applied to the study of cosmic neu¬ 
trinos. The consistency of vacuum oscillations and IceCube observations was discussed in 
m, m, |21] compare with [22]. In [20] the impact on the flavor ratio of the uncertainties 
on oscillation parameters was analyzed and the present work develops the discussion. 

In this work, we show that the parameterization of neutrino oscillations in vacuum intro¬ 
duced in m can be promoted to an exact parameterization and can be argued to be the 
natural parameterization for the discussion of oscillations of cosmic neutrinos)]] We obtain 
new and exact expressions of the three parameters in terms of the known mixing angles and 
CP violating phase. We evaluate their numerical values and uncertainties, illustrating their 
usefulness by discussing three applications: 1) we compare the predicted flavor fractions 
and those that are allowed by the present observations; 2) we quantify the uncertainties in 
the prediction of the fraction of muon neutrinos due to oscillations; 3) we argue that, even 
after accounting very conservatively for the uncertainties of oscillations, the intensity of 
the Glashow resonance |24] differs greatly in the alternative cases of pp- and py-production 
as remarked in [TO] and later discussed in mmmm- Throughout this work, we argue 
that a Gaussian treatment of the errors of these natural parameters is quite adequate for 
the present precision. 


2 Natural parameters for three flavor vacuum oscillations 

In this section, we motivate, define and analyze a parameterization of vacuum neutrino 
oscillations, elucidating the relationship between this and other parameterizations. 

We begin by counting the number of independent vacuum oscillation parameters [T] 

n 

Pit' = \ U ti\\ U i'i\ where t = t... (1) 

i= 1 

in the case of n light neutrinos. 

The vacuum oscillation formula depends upon the squares of the leptonic mixing matrix 
\UJil. They correspond to (n — l) 1 2 independent parameters, as it is clear considering all 
| Ugj | as independent (albeit constrained) parameters except the ones of the first row and 
column, that can be obtained from unitarity, e.g., \U^\ = 1 — X^=2 Wti\- 

1 An alternative choice of the parameterization that gives different insight on the allowed ranges of the 

oscillation probabilities is discussed in m- 


2 



But the vacuum oscillation probabilities P'w are also symmetric in the exchange of the 
flavor indices £ Pr £!, thus they require less parameters. Since we can again rely on similar 
unitarity relations when we sum on all flavors, namely Pp^ = 1, this implies that the 
number of independent parameter is just n(n — l)/2. This means that when n = 3 we have 
3 independent parameters, when n = 4 we have 6 of them, etc. 

From here on and in view of the present experimental situation, we focus on the three 
flavor case (n = 3) where we have three parameters, as first remarked (to the best of our 
knowledge) in ref. [HE] . 


2.1 Motivations and definition 


Next, we motivate and introduce the choice of the three natural parameters. The param¬ 
eters Pq, P \, P 2 are defined as follow, 


P — 4 

D _ ree 3 
P) — -r- 


, Pi = 


P — P 

1 e/i 1 er 


, P2 = 


P 1141 + Ptt — 2P L 


flT 


( 2 ) 


We can write in terms of Pq, Pi, P2 the matrix that contains the probabilities of oscillations 
of cosmic neutrinos. This is the following symmetric matrix, 

/ 1 + 2 Po 1 -Pq + P 1 l-Po-Pi \ 

V=\ 1 + f - P, + P 2 \ + ^~P2 (3) 

V | + f + Pi + p 2 / 


It acts on the vector of fluxes before oscillations P° = (Fg, F^, F°) just as F = V F°, 
giving the vector of fluxes observed after oscillations, F = (F e , F M , P T ). 


We would like to give our reasons (that are largely based on the available experimental 
information) why we consider that Eq. [ 2 ] is the optimal choice of parameters. 

1. The oscillation probability that is singled out in Po is P ee , that is well-known (be¬ 
ing directly measured by low energy solar neutrino experiments and probed also by 
reactor and high energy solar neutrino experiments). 


2. The difference of P efl — P eT contains most of the uncertainties. 

3. The last combination of oscillations probabilities, P 2 , is positive and very small. 

4. A specific choice of the overall coefficients is adopted in order to have coefficients that 
are either zero or close to 1 in the expressions of all oscillation probabilities, Eq. [3] 

5. Setting Po = Pi = P 2 = 0, all oscillation probabilities become Pw = 1/3, namely, 
any information on the original flavor is lost: The three parameters describe the po¬ 
tentially measurable information on flavor that survives cosmic neutrino oscillations. 
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Parameter 

Mean value 

Standard deviation 

Po 

0.109 

0.005 

Pi 

0.000 

0.029 

P 2 

0.010 

0.007 


Table 1: Table of present values and errors of the natural parameters. 


The above argument clarifies that it is possible to introduce such a parameterization di¬ 
rectly, without the need to associate it to a specific scheme of approximation of the mixing 
matrix, but rather, keeping it exact. However, various approximation schemes used in 
the literature have allowed to uncover the most interesting properties: The second one is 
known since hd while the first evidences of the third one were found in m- A detailed 


comparison with other parameterizations used in the literature is offered in Sect 2.4 


2.2 Connection with the standard parameters of neutrino mixing 


Compact and useful expressions of the natural parameters in terms of four standard param¬ 
eters, the mixing angles @i 2 ,# 23 ,#i 3 and CP phase violation phase 5 , are as follows, 

sin 2 6*12 ’ 




1 - 


2 1 

+ £ "3 


1 - e 


7 cos 26*12 + P 


1 -3e 


Pi = 


P2 = 2\ 72 + 4 ^ (1_e) 
where we introduce for convenience the following 4 small parameters, 
e = sin 2 @13 a = sin @13 cos S sin 26*12 sin 2@23 

P = cos 26*23 7 = a — f cos 2 @i 2 (l + e) 


(4) 


(5) 


These expressions are new and exact. Note property 3 listed in Sect. 2.1 of this parame¬ 
terization. 


2.3 Numerical analysis 

The parameters a, P, 7 are small and to date not known precisely, whereas e is very small 
and precisely known. We can then order these parameters according to their (presumed) 
size, and consider sin 2 26*12 as zeroth order; sin @ 13 , cos 26 * 23 , oc, P, 7 of first order; e of second 
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Figure 1: Distribution of the natural parameters Pq, Pi and P 2 , due to the uncertainties in the 
mixing angles and the phase of leptonic CP violation. 


order. In the same sense, Pq is a zeroth-order parameter; Pi is first order in a and /3; P 2 
is second order in a and (3. Note that P 2 is bound to be positive. 

Using the present knowledge of mixing angle and CP violating phase [25], we obtain the 
values and the errors of the natural parameters. We show the results in the Table [lj 
assuming normal mass hierarchy. It is easy to repeat the same steps with inverted hierarchy, 
but the differences are not large. From this Table we notice that with present data the 
average values obey (Pq) (Pi) ~ (Po) whereas their variances obey 5P\ 5Pq — 5 P2- Pq 
is well known, because is related to survival probability of solar low energy neutrinos and 
#13 or e is well measured by reactor experiments. As we see from figure [lj Po and Pi are well 
represented by Gaussian functions; P 2 is not Gaussian but it is a very small parameter. For 
these reasons, as we argue in the rest of this work, we can use a Gaussian approximation 
without introducing severe inaccuracies in the numerical analysis of the oscillations. This 
is a new result, that allows one to obtain convenient analytical expressions for different 
examined quantities and to quantify easily the uncertainties. 

The probabilities of oscillation given in Eq.[3]have a very simple form: they depend linearly 
upon the natural parameters. Moreover, in first approximation, they could be expressed 
only in terms of Po, because Pi and P 2 give small corrections. Using the value of Table 
[I] and the natural parameterization of oscillation matrix, we obtain the probabilities of 
oscillations, 

Pee = 0.552 ± 0 . 010 , Pep = Per = 0.224 ± 0.029 . . 

P[it = 0.378 ± 0.008, P /t/i = P TT = 0.398 ± 0.029 ^ ’ 

Two couples of probabilities have (almost) the same values, because with the present best 
fit value (Pi) = 0 and the numerical differences between these expressions are small. 
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2.4 Comparison with other parameterizations 


First of all, we consider the leading order in the small parameters $13 and cos 2#23- The 
parameter y that was introduced in HD, namely 


1 


1 


y = sin 20 i 2 cos 2023 H— sin 2012 cos 2012 0 i 3 cos 5 


(7) 


coincides at this order with Pi, while higher order terms are neglected. In HD, the errors 
on P eM and P MM were estimated to be 0.05; after 10 years, these errors amount to 0.03 and 
we see from table [l] that, still, the uncertainty is mostly due to Pi. Equivalent parameters 
have been used also by other authors: E.g., in ESI, we read that 


A “universal” parameter related to B has been noted previously in the literature. 
It is called —A in Z. Z. Xing, Phys. Rev. D 74, 013009 (2006), and +A in W. 
Rodejohann, JCAP 0701, 029 (2007). 


although, curiously, reference HD is not mentioned. Note that the parameter y or equiva¬ 


lently A satisfies properties 2 and 4 of Sect. 2.1 


As already mentioned, a three-parameters description of the probabilities of vacuum oscil¬ 
lations with was first introduced in ESI, where the parameters A, B, C were defined: see 
Eq. (19) and note the symbol 0(e 3 ) used there to emphasize the use of quadratic expression 
in the small parameters that quantify the deviation from tribimaximal mixing, collectively 
denoted as e. This approximation is even better than the linear approximation, and thus 
was argued to be adequate for the present needs m- However, there are two evident short¬ 
comings in the type of procedure adopted there to introduce the new parameterization: 
The emphasis on tribimaximal mixing given in Eq. (7) and see Eq. (16) of reference |16| is 
felt as artificial to date, especially now that the measurements showed that 0 i 3 is non-zero, 
contradicting the most interesting prediction of tribimaximal mixing. Moreover, the efforts 
used to obtain a quadratic the expansion is also unnecessary: the parameterization can 
be promoted without significant efforts to an exact one, as the expressions in Sect. | 2 . 2 | are 
valid to all orders in e and easy to use. Thus, a direct, valid and advantageous procedure 
is to introduce the parameterization since the start, as in Sect, 
with our parameterization is, 


2.1 The detailed relation 


A = 9P 0 -1, B = 18 Pi and C = 18 P 2 . 


( 8 ) 


Both of them share features 1, 2, 3 discussed in Sect. 2.1 and they can be termed as natural. 
In the following we adopt Pq,P\, P 2 due to features 4 and 5 of Sect. |2.1| and to the fact 


that symbols reflect the hierarchy noted in Sect. 2.3 The parameterizations are however 


equivalent and it is easy to compare the results obtained with them. 


Another equivalent parameterization was used in m This begins from the “universal” pa¬ 
rameter of HD and improves the description of the oscillation probabilities by introducing 
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a new parameter. Also this parameterization is introduced in connection to tribimaximal 
mixing and using the same quadratic expansion of the previous parameterization: note the 
symbol ~ of approximate equality in Eq. (2.5) of |17| . Thus, the same comments on the 
methodology made just above apply also to this case. The relation with our parameteri¬ 
zation is, 

A = Pi and ^- = P 2 (9) 

A third parameter is not introduced, being replaced by 0\ 2 and \U e s\ = sin$ 13 . Properties 

2 , 3 and 4 of Sect. |2.l| are all satisfied. Therefore, assumptions and results can be compared 

easily: e.g., Tab. [I] implies P 2 < 0.017 (resp., 0.031) at 1 sigma (resp., 3 sigma) whereas 

the value given in Eq. 16 of E! implies P 2 < 0.0145 (resp., 0.0465). Similarly, the values 
—2 

A ~ 0.02 and A ~ 0.008 quoted there correspond to P\ — 0.02 and P 2 ~ 0.004, that 
are included in the 1 sigma range of Tab. [lj The differences are due to the improved 
measurements of the oscillation parameters since 2007 and in particular to the inclusion of 
#13 that now is measured and known to be non-zeroj^] Note finally that in Fig. 1 of [[16] 
the expected ranges of the parameters are presented, and these can be compared with our 
Fig. [I] and Tab. |T} 


In the literature also other different linear combinations of the parameters have been con¬ 
sidered, see e.g. [23]. The parameters X,Y,Z introduced there, however, do not satisfy 
features 1, 2 , 3 and 5 of Sect. 24, and in this technical and restricted sense, we do not 
call them ‘natural’. The choice of parameters in [23] has its own motivations but it is less 


useful to keep under control the impact of the uncertainties on oscillations that is one of 
the main goal of the rest of this work. 


3 Applications 

We will consider two quantities, that are affected by oscillations; the flavor ratios and the 
fraction of events due to Glashow resonance. 

We denote the fractions of i/£ at source and the one at Earth (i.e., after oscillations) 
respectively as, 

0° = Ff/ E F e and & = F tt E F * ( 10 ) 

t a 

where of course Yle F e = Yle F t- Suppose that the initial flavor ratio is given by 

(&&£) = (! ~9-h,g,h). (11) 

2 Let us repeat that this differs from what was expected from the tribimaximal mixing scheme, that has 
been emphasized in nni and na. 
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Figure 2: Flavor triangle. The present observational information on the flavor composition of 
cosmic neutrinos at 1, 2 and 3a analyzed assuming a = 2 is compared with the expectations derived 
for four different hypotheses on the mechanism of production of the neutrinos. The ellipses derive 
from a Gaussian treatment of the error based on Ta&ZeQJ see also Appendix [~T| 


After propagation the flavor ratio is modified as follows, 

= 5 + ( 2 ~ 3 g - 3h)P 0 + (g~ h)P\ 

| + ^( — 2 + 3 g + 3h)Po + (1 — 2 g — h)P\ + (g — hfP^ (12) 

£r = g + ^( — 2 + 3 g + 3/r)-Po + (~ 1 + 9 + 2 h)Pi — (g — h)p 2 

Below, we will emphasize £ /t since it is quite directly connected to an observable quantity, 
namely, the fraction of track-type events. 


3.1 Flavor ratio after the oscillations 

A first application of the natural parameterization is the study of the flavor ratio of neu¬ 
trinos, considering different mechanisms of production]^] The impact of uncertainties on 
the mixing angles and CP violating phase on the triangle of the flavors was first discussed 
in Here we will update the analysis by using updated values of the uncertainties. We 
will verify that a Gaussian treatment of the natural parameters leads to results in good 


3 Intermediate possibilities have been also considered in |25l .29 . 
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fraction of e : fraction of mu 


pions decay 
charm 
neutrons 
damped 


Figure 3: Same as previous figure but using the spectral index a = 2.3 in the analysis of the IceCube 
data. 


agreement with Monte Carlo simulation and, moreover, we will compare our theoretical 
analysis with three years IceCube HESE data + passing muons. 

We consider: 

1. pions decay (h = 0 ,g = 2/3; blue); 

2 . neutrons decay (h = 0 ,g = 0; green); 

3. damped muons {h = 0,g = 1; red); 

4. charm mesons (h = 0, g = 1/2; orange). 

Using Eqs. ( |12[ ), we represent the allowed regions by propagating the errors on the predic¬ 
tions by a Monte Carlo simulation; this gives the 4 dotted regions of the flavor triangles in 
Fig. [2] and [3| 

These regions can be compared with those obtained with a Gaussian treatment of the 
errors on Pq, Pi and P 2 . Following the implementation of Appendix [Aj we obtain the 
4 ellipses of Fig. [ 2 ] and [5] that enclose the 99% CL regions. We see that the differences 
between these the Montecarlo and the Gaussian treatments are not very important. The 
Gaussian approach seems to be appropriate for the present needs. Note that the latter is 
significantly easier to implement. 

Let us repeat that the expected theoretical regions (the dotted areas obtained by Mon¬ 
tecarlo and the elliptic curves corresponding to Gaussian approximation) depend only by 
initial flavor ratio and are not affected by the energy spectrum of the neutrinos, that is as- 
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sumed to be universal for all neutrinos. In other words, the four theoretical regions shown 
on the flavor triangles of Fig. [2] are just identical to those of Fig. [3j 

On the contrary, the 1 <7, 2 a and 3er zones do depend on the energy distribution of the 
neutrinos. The confidence levels, indicated in the flavor triangle, correspond to the result 
of the data analysis of IceCube events (the high energy starting events, whose initial vertex 
is in the detectors and the passing muons, i.e. the thrugoing muons) discussed in |21] and 
They have been obtained assuming a power law distribution given by, 

d(j> 


dE„ 


— <i>oE v ' 


(13) 


Thus, the observed flavor ratio is a function of spectral index a. In Fig. [2] we have used 
the value preferred by the simplest theoretical expectations, namely a = 2. 

For comparison, we have shown also the case a = 2.3 in Fig. [3j namely the best fit value 
of the dataset of high energy starting events collected by IceCube in the first three years. 
The comparison with Fig. [2] shows that the actual value of the spectral index plays some 
role in determining the allowed regions: the steeper spectrum a = 2.3 enhances the role of 
electron neutrinos and diminishes the one of muon neutrinos (whose effective area is very 
small close to the threshold) thus it requires to increase the content of muon neutrinos at 
the source in order to reproduce the observed track-to-shower ratio. For this reason, 
the agreement of the neutron decay scenario with the data worsens for a = 2.3. However, 
this kind of effects it is not yet crucial for the analysis and in particular the neutron decay 
scenario is not yet excluded. In fact the most important conclusion is just that the small 
number of events presently available, does not allow us yet to exclude any mechanism of 
production [20]. This remains true also using a = 2.6, namely the best fit including also 
low energy events. 


The confidence levels are in reasonable accordance with those of IceCube data analysis 
(see again figure [2]) and with those of [20] . The uncertainties due to the oscillation pa¬ 
rameters have been presented in m in a different manner, but the results are in excellent 
agreement^] 


3.2 Errors on flavor ratio 


Let us discuss further the point of the errors. The flavor ratios given in Eqs. [12] depend 
linearly by Pq, P\ and Pi. Therefore, it is straightforward and quite easy to evaluate the 
Gaussian errors on these quantities. Let us focus on £ /t . 


4 A recent work ;32] appeared after the present one also shows the uncertainties due to oscillations in 
the flavor triangle. The Monte Carlo procedure is used and the results coincide with our ones. 
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Figure 4: Gaussian errors on the fraction of on Earth at 1 sigma, as a function of the neutrino 
fractions at the source. 

From the formula for ^. again Eq. ( |12[ ), we see that the term linear in P\ becomes very 
small when v e and are about equal (e.g. charm mesons), and Fig. [4] confirms that that 
this type of mechanisms gives very small errors on the flavor ratio measured at Earth. 
Indeed, an initial composition of 1:1:1 would not be modified, or in other words, the error 
would be just null). On the contrary, the mechanisms that produce only z/ e (neutrons decay) 
or only (damped muons) give the biggest error, about 10% on final flavor ratio. The 
pion decay, that is the most plausible mechanism, is between the two extreme situations; 
the error on the muon fraction is about 3%. 

Let us remark that despite the relatively large uncertainties on 023 and 5 oscillations, 
the uncertainties on % are small also in the worse scenario, namely, the neutron decay 
scenario. 

These results are in good agreement with those of m, see in particular Fig. 4 there. 
Moreover, we note that the expectations from the pion decay mechanism agree quite well 
with the results of the analysis of the existing data; see again Fig. [2] 


3.3 Glashow resonance 

With the formalism of this paper it is easy to write analytical expressions for some inter¬ 
esting signal including the effect of three flavor oscillations. We analyze the case Glashow 
resonance [24], i.e. the production of W~ starting from electron antineutrino due to the 
process, 

V e + e~ —> W~ (14) 
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Figure 5: Fraction of v e due to pp or py interaction. The distribution obtained by a Montecarlo 
extraction compares well with the Gaussian distributions obtained from Table^ (continuous lines). 


The process is possible when the antineutrino with energy greater than 6.3 PeV collides 
with an electron at rest. 


An interesting point for us is that different astrophysical mechanisms produce different 
fraction of z7 e , as already noted in [IQj and further discussed in mmmm- From 
the pp interaction, e.g. ED], we can obtain all the type of pions, instead from py interac¬ 
tion, e.g. [31], we obtain mostly 7r + and 7r°. After the decays the flavor ratio for pp are 
approximatively equal for neutrinos and antineutrinos, 

£ = (l,2,0)/6, £° = (l,2,0)/6 (15) 

while for py, the neutrino and antineutrino channels contain a different number of particles 
and lead to different flavors, 

£° = (l,l,0)/3, £g=(0,l,0)/3 (16) 


where we have normalized the two fluxes to a single particle. In the case of py interaction 
v e are not produced at the source, but only after the oscillations. The fraction of electron 
antineutrinos at Earth are given by a linear expression in the parameters Pq and Pi- 


cPP 

s>u e 


cPl 

Sy e 


1 

6 


+ 


1 

3 


Pi 


1 — 3Pq 
9 


+ l Pl 


(17) 

(18) 


These two distribution can be obtained with Monte Carlo simulation, using the distribu¬ 
tions of mixing angle and CP phase violation. At 3a (i.e., 99.7%) we find that 


tvp 

cP 7 

St/g 


0.167 ±0.029 


0.075 ± 0.029 
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where the uncertainty is due to Pi. These ranges compare well with the maximum and 
minimum values given in Fig. 6 of m- Note that we are using the newest measurements 
of the oscillation parameters and in particular we have included the effect of $13, that is 
now measured and known to be non-zero. 

If the total flux of neutrinos is the same in both cases, this implies that the signal from 
the Glashow resonance is about two times weaker in the py case than in the pp case. 
Thus, oscillations give a final flavor ratio that is significantly different in the case of py 
production mechanism. The results are further illustrated in figure [5} Note that, also for 
this application, the very small differences with the result of the Gaussian approximation, 
evident from Fig. [5j justify the use of a linear approach. 

If we compare the two cases assuming to have the same cosmic ray flux and the same 
amount of collisions, as done in the difference between the two cases is even more 
dramatic. In fact, the py mechanism leads to neutrinos of lower energies than those from 
the pp mechanism, and thus the chances of observing a neutrino signal from the Glashow 
resonance decreases. 


4 Summary and discussion 


The results of IceCube have greatly increased the interest in an accurate description of 
propagation of cosmic neutrinos, accounting in particular for the minimal hypothesis of 
three flavor oscillation in vacuum. In this work we have discussed what is the natural 
choice of the parameters to describe vacuum oscillations of cosmic neutrinos. 


We have shown that the parameterization of neutrino oscillations in vacuum introduced 
in [il6] can be promoted to an exact one (without any need of performing expansions, 
approximations or of making references to case of tribimaximal neutrino mixing) and can 
be argued to be the natural parameterization for the discussion of oscillations of cosmic 
neutrinos, see Sect. 2.1 We have obtained new and exact expressions for the three relevant 


parameters in terms of the standard mixing parameters, see Sect. 2.2 


In Sect. [3j we have illustrated the usefulness of the natural parameters Po, Pi, P 2 given in 
Eq. [2] by discussing the expectations on the neutrino flavor ratios and their errors. We 
have also analyzed the expectations on the intensity of the signal due Glashow resonance, 
that depends on the mechanism of neutrino production. We have included the effect of 
the uncertainties on oscillation in the flavor triangle, comparing the predictions with the 
results of the analysis of the flavor of cosmic neutrinos seen by IceCube, finding results in 
good agreement with the previous literature, [20] and [21]. We have confirmed Em that the 
pp and py production mechanism lead to significantly different predictions for the intensity 
of the Glashow resonance. 
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With improved data and analyses of IceCube data, it will be more and more important to 
include the effect of the various theoretical uncertainties, including those due to parame¬ 
ters of neutrino oscillations. Here, we have shown that the expectations obtained with a 
Gaussian treatment of the natural parameters Po, Pi, P 2 are very similar to those obtained 
in more complete descriptions of three flavor oscillations baseed on |20| . We have proved 
that the Gaussian treatment, particularly easy to implement, is quite adequate for the 
present needs. 

Let us conclude by stressing that the parameters Po, Pi and P 2 can be used also to provide 
us with an ideally exact and compact description of cosmic neutrino oscillations. E.g., the 
distributions of these parameters obtained with MonteCarlo methods, described above and 
illustrated in the flavor triangles of Figs. [2]and[3j go beyond the Gaussian approximation 
and it is formally accurate. 

Future global analyses of the data on three flavor neutrino oscillation will be able to 
derive which are the precise distributions of the natural parameters of the cosmic neutrino 
oscillations and their correlations, simply because these parameters are functions of the 
conventional mixing angles and CP-violating phase. In view of the above discussion, it will 
be particularly useful to have precise distributions of the parameters Pi and P 2 . 
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A Allowed regions in the Gaussian approximation 


Let us consider the two dimensional Gaussian likelihood, 


c = exp \—\{v — (v}Y S 2 (v — (11))] 

27r-^/det(E 2 ) 

A confidence level (0 < C.L. < 1) defines the allowed region C > £ max (1 — C.L.). Its 
contour is an ellipse that can be obtained from the following parametric expression, 


where v = 


E 2 = 


u 2 , 2 


( 21 ) 


\y(p)J 


(^ ) ) )+ v /-21og(l-C'.L.) 


cos 9 sin 9 
— sin 9 cos 9 


a + cos (p 
O- sin</? 


( 22 ) 
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where </? = [0, 2 i r] and where we defined, 

„ 1 f 2cr 2 


6 = - arctan 
2 




(23) 



a x + °y i (o'y ~ crDVl + tan 2 26 


In the flavor triangle, we have known linear combinations of Pq,Pi,P 2 , 


x = (f M - £e)/ \/3 = x 0 + XiPi and y = £ T = y 0 + ) 

From Tab. [lj one evaluates (x) = xq + Xi(Pi), a 2 = x\ SPf, o 2 = x^y* 5P 2 , etc. 


(24) 
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